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The simple microscopic model for the high-Tc superconductors based on the electron-phonon 
(EPH) and electron-electron-phonon (EEPH) interactions has been presented. On the fold mean- 
field level, it has been shown, that the obtained model supplements the predictions based on the 
BCS van Hove scenario. In particular: (i) For strong EEPH coupling and T < Tc the energy gap 
(Atot) is very weak temperature dependent; up to the critical temperature Atot extends into the 
anomalous normal state to the Nernst temperature, (ii) The model explains well the experimental 
dependence of the ratio Ri = 2A(JJj/fcsTc on doping for the reported superconductors in the terms 
of the few fundamental parameters. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.20.Mn, 74.25. Bt, 74.72.-h 



I. INTRODUCTION 

In the study, we present the simple microscopic theory of the high-Tc superconductivity [l[ . The organization of 
the paper is as follows: 

In Section II we call for the pairing mechanism. First of all, we discuss the main experimental and theoretical 
results. Next, on the basis of the presented analysis, we give the postulates, which determine the microscopic model 
for the high-Tc superconductors in the second quantization form. 

In Section III, by using the unitary transformation, we deduce from the initial Hamiltonian the simple mean-field 
model. In the framework of the toy model (only the s-wave state) , we discuss the properties of the energy gap in the 
superconducting and Nernst region; the numerical predictions are supplemented by the analytical approach. Next, for 
selected high-Tc superconductors, we calculate the fundamental parameters of the model and compare the obtained 
theoretical results with the experimental data. 

Finally, in Section IV we summarize the results. 

The main body of the paper is supplemented by the Appendixes. The formally exact expression for the self-energy 
matrix is obtained in Appendix A. In Appendix B we present in detail the fold mean-field approximation. The van 
Hove and generalized mean-field thermodynamic potentials are calculated in Appendix C. The lists of the experimental 
values of the thermodynamic parameters for the selected high-Tc superconductors are collected in Appendix D. 

II. THE PAIRING MECHANISM 

The real cuprates are three-dimensional. However, their physical properties are strongly anisotropic. On the basis 
of very small coherence length in the c direction (smaller than the interplane distance) one can suppose that the 
Cu02 electrons play the special role in the physics of the high-Tc superconductors [11 . Unfortunately, the pairing 
mechanism for the planar problem remains highly controversial and many different hypotheses are suggested. In the 
literature two fundamental directions in search for the pairing mechanism have been crystallized. The first approach 
is based on the single-band Hubbard model, its extensions or related models e.g. t — J model the second 

approach emphasizes the relevance of the electron-phonon interaction ^Qj. 

Why is the Hubbard model so much studied? First of all, some analysis suggest that the one-band Hubbard model 
reproduces well the spectra of the more complicated three-band Hamiltonian for electrons in the copper oxide planes 
(the Emery model) Q, Q. For example, by using the finite cluster method, Hybertsen et al. have shown that 
the one-band Hubbard model with the small next-nearest-neighbor integral t should have the following parameters: 
t = 430 meV, t = — 70 meV and Uh = 5.4 eV [lO|. Secondly, for the half-filled electron band and large on-site 
Coulomb interaction, the Hubbard model reduces to the Heisenbeg model, which describes well the spin dynamics of 
the underdoped high-Tc superconductors [2]. On the basis of the quoted facts some authors suppose, that the strong 
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FIG. 1: The graphical representation of the electron-electron-phonon coupling. The electron states are represented by the 
straight lines, while the phonon state by the curly line. The vertex is denoted by the black dot and its strength is given by 
(?'^' , (q, 1). In Fig. [T](A) the electrons emit the phonon during the scattering; in Fig. [T](B) the electrons absorb the phonon. 



electronic correlations modeled by the Hubbard model can alone induce the superconducting state in the cuprates. 
Unfortunately, the studies carried out by several groups have shown that the Hubbard model gives no obvious evidence 
for superconductivity with the large critical temperature [111] . On the other hand, there is the strong tendency for 
superconductivity in the attractive Hubbard model for the same value of the on-site Coulomb interaction. Finally, we 
noticed that probably also the three-band Hubbard model and the t ^ J model do not superconduct at temperatures 
characteristic for the cuprates (l^ . 

The relevance of phonons to the pairing mechanism in the high-Tp superconductivity also constitutes a complicated 
problem. On the one hand there exist many experimental observations which have been taken as evidence for the 
electron-phonon interaction in the cuprates. For example: the strong isotope effects on Tc in the underdoped 
superconductors p^ . the phonon renormalization in the Raman measurements [l^ . the phonon- related features of 
I — V characteristics obtained by using the tunnelling experiments [l^ and the dependence of the penetration depth 
on the substitution O^^ by O^* [lH. Especially important results come from the ARPES measurements which give 
the evidence on the low-energy kink in the quasiparticle spectrum around the phonon energy both for the nodal and 
antinodal points [13], ^M'l ^^^o the ARPES isotope effect in Re (S) has been observed [19[. On the other hand the first 
principles calculations support the view that the conventional electron-phonon coupling is small [lO]. For example, 
Bohnen et al. have predicted that the electron-phonon coupling constant for YBa2Cu307_y is equal to 0.27; so the 
strong Hubbard correlations should completely suppress the phonon- mediated superconductivity |2l| . 

After summarizing the mentioned experimental and theoretical results, one can conclude that: (i) the cuprates 
belong to the strongly correlated systems but probably these correlations in the superconductivity domain are beyond 
the Hubbard or related approaches, since in these models the pairing correlations are too small, (ii) in the cuprates 
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the conventional electron-phonon interaction is small but according to the experimental data one can suppose that 
the phonons play the important role in the pairing mechanism. 

In order to solve the problem of high temperature superconductivity we present and examine the following scenario: 

(i) In the superconductivity domain of the cuprates the fundamental role is played by the electrons on the Cu02 planes. 

(ii) In the cuprates exists the conventional electron-phonon interaction, which has not to be strong. 

(iii) In the cuprates exist strong electronic correlations, but the electron- electron scattering in the superconductivity 
domain is inseparably connected with the absorption or emission of vibrational quanta. 

In the simplest case the first and second postulate coincides with the phonon-van-Hove-scenario for high-Tc su- 
perconductors [131, [l^. The third postulate states that the effective electronic correlations in the superconductivity 
domain are connected with the three-body process: the electron-electron-phonon interaction. In Fig. [T] we show 
in detail the diagrammatic representation of this interaction. We notice that the EEPH coupling has a significant 
property which distinguish it from the Hubbard interaction; it does not destroy the classical phonon-mediated pairing 
correlations. Additionally, one should pay attention to the fact, that the first postulate has also the essential signifi- 
cance for the third postulate. Namely, for the two-dimensional case (the van Hove singularity at the Fermi level), the 
EEPH coupling has significantly strong influence on the physical properties of the system (see next section). 

Below, we consider the Hamiltonian that describes the postulated coupling of the correlated electrons to phonons 
in the second quantization form: 

i? = i/(")+H(i)+i7(2). (1) 
The first term represents the non-interacting electrons and phonons: 

H(o)^^ekcLck. + 5]a.q6t5^, (2) 

ka q 

where Ek = £k ^ ^k and /i denotes the electron band energy and the chemical potential respectively. For the 
two-dimensional square lattice and the nearest-neighbor hopping integral t, we have: £k = ~^7(k), where 7 (k) = 
2 [cos (kx) + cos (ky)]. The symbol Uq stands for the energy of phonons. The interaction terms are given by: 

^''^=Effk^(q)4+q.Ck.0q, (3) 
kq<T 

and 

H''' - E C' 1) 4-i.Ck.ct,^,^^_^c,,_^^, (4) 

kk qlcr 

where 0q = feLq + bq. The matrix elements g^"^ (q) describe the electron-phonon interaction (2^ and the symbol 
g^ ^, (q, 1) determines the strength of the electron-electron-phonon coupling. Since the Hamiltonian is the hermitian 
operator, we have: 5*'^^ (q) = g^^l^ (-q) and g^^^, (q, 1) = S^'+i+q^k-i ^ + '^)- 

III. THE FOLD MEAN-FIELD THEORY 

From the mathematical point of view, the Hamiltonian ([!]) belongs to the class of so-called dynamic Hubbard models. 
These Hamiltonians describe modulation of the Hubbard on-site repulsion by the boson degree of freedom and were 
studied in detail by J.E. Hirch, F. Marsiglio and R. Teshima in the papers p5j-[27|. 

In the simplest case, the dynamic Hubbard Hamiltonian describes the Holstein-like electron-boson interaction, where 
the coupling constant increases with the electron occupation. This model, however interesting, is inappropriate in the 
description of the high-Tc superconductivity. In the more elaborate scheme, the properties of the dynamic Hubbard 
model, can be studied by using the generalized Lang-Firsov transformation: H — e^ IIe~^ , where the generator S 

has the form: S = g (^bj — (rij-f- + nji — Uj-^nj^) and Uj^r is the number operator for the site j and the spin a. 

Next, the obtained effective Hamiltonian should be analyzed in the framework of the Eliashberg-like approximation 
[131- [2^; the central result is that the critical temperature increases with retardation. 

We notice that the form of the Hamiltonan ([1]) suggests using the Eliashberg approach in order to analyse the physical 
properties of the considered system [29']. This method is appropriate, because one can retain simultaneously the 
electron and phonon degrees of freedom. Additionally, in the framework of the Eliashberg formalism, the formally exact 
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FIG. 2: The four-body scattering event contributing to the interaction part of the Hamiltonian H10[). The oval connecting the 
electron's lines is an illustration of [/, , / (q, 1) which denotes the effective potential. 



expression for the self-energy matrix is easy to determine. The analysis of the above issue is presented in Appendix 
A, where we have shown that: (i) in the case of the absence of the lattice distortion, the first-order contribution 
to the self-energy is equal to zero, (ii) the electron-phonon coupling is directly increased by the electron-electron- 
phonon interaction, (iii) the form of the second-order contributions to the matrix self-energy is very complicated. 
Unfortunately, this fact hampers considerably the detailed calculations on the Eliashberg level. 

In the case of the presented paper, we have eliminated the phonon operators only from the EEFH term. The 
obtained Hamiltonian is next added to the BCS operator [s^. As discussed in the following subsections, presented 
approximate description represents a generalization of the BCS van Hove scenario [23l |. 



A. The unitary transformation 



In order to eliminate the phonon degrees of freedom in the Hamiltonian ^ we use the Frohlich-type unitary 
transformation: 



where the operator S denotes the generator: 



i?' =e-*^ (ij(0) + ^(2) j ,.5 



5 _ ^ S'q 



(5) 
(6) 



In our case: Sq = 7qfeq -I- 'ylbt, where 



7q= XI * (k,k',l,q) 4_,^Ck<xc|^ 



k +l+q-(T 



kk la 



The expression ([5|) can be rewritten in the approximate form: 



(7) 



(8) 



+ i 



where the square brackets [, ]_ denote the commutator. To eliminate the second term in Eq. (|S]) we assume that the 
generator fulfills the relation: i[H^^\S]- + H'-^^ = 0, hence: 



$(k,k',l,q) = 



^^il (q,l) 



£k-i - ffk + '^k'+i+q - ^k' - '^q 



(9) 



Next, we can reduce the Hamiltonian H to the following expression: 
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Heff =< Ophl^^' |Oph >^ ^kcLck^T + ^k,k' 1) cLi<TCk<Tcj^'^j^q_^ v_^ct_i^,_,_^^c_k'^cLk+i_^c_k-^, 

k<7 kk'qlff 

(10) 

where the phonon vacuum state is given by |Oph >• The pairing potential (q, 1) is of the form: 

t/,,.(q,l).- (11) 

(^ek-ek-i + ^k' -Ek'+i+qj -^0 

We assume additionally that: g^^^, (q, 1) ~ 5'^-* and Uq ~ wq; the symbol ujq denotes the characteristic phonon 



frequency. On the basis of the expression pOj) we conclude that the EEPH interaction can be replaced by the effective 
four electron-electron (4EE) scattering event; the diagrammatic representation of this interaction is shown in Fig. [21 
From the Eq. ([TT|) it is clear that the effective potential is attractive if: |ek — Ck-i + Ck' ~ ^k'+i+ql '^o- In the 
subsection, we consider the simplest case, when the attractive part of the 4EE potential can be written as: 

for |£k| < 1^0 • We have used the factor because the potential energy term represents the interaction between every 
four of particles counted once. 

By using the fold mean- field (MF) approximation the Hamiltonian ^TU\\ takes the form (see also Appendix B): 

Heff - E^'^^.Ck. (13) 

uio 

■ Y I A.P (Aact.c^-k. + A;c_k.ck-.) . 

kfT 

The symbol J^k" denotes the sum over the states when the 4EE potential is attractive; A^- = X^k" (c-ko-Ck-o-)- 



ko- 

Ulo 



u 
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B. The toy model 

The thermodynamic parameters of the high-Tc superconductors can have essentially different properties in compar- 
ison with the low- Tc" materials. From this reason, the analysis of the results obtained in the framework of the simplest 
approach is important, because these predictions facilitate the interpretation of the fundamental experimental data. 

Taking into account the operator (jl3p we can write the total Hamiltonian in the form: 

H''"^ ^Y'kcl^'^. (14) 

ko- 

^ (v+ ^\aA g (Act,cL,, + A-^c-kiCkt) , 

k 

where we have omitted the chemical potential and A = A^. In Eq. 1^^, the symbol V represents the BCS pairing 
potential obtained from the Hamiltonian Q. 

Now, we establish the energy scales in the presented model. The nearest-neighbor hopping integral t is of the order 
of (200 — 400) meV [3l|-[33|- We notice that in the numerical model calculations we take t as an energy unit. From 
the ab initio calculations arises the fact that V £ (0, 2t) |2l|. In order to determine qualitatively the possible values 
of U, we note, that the simple BCS pairing potential is obtained from the expression: 

-I"'"!' ^ " (15) 



(fk-fi+,)'-"l 2iV 

where the characteristic phonon frequency ujq is of the order of Debye frequency {uju ^ 0.3t); the electron-phonon 
coefficient g^^^^ has nearly the same value [35|. Next, we assume that the largest energy in the high-Tc superconductors, 
of order (5 — 10) eV, is the electron-electron-phonon potential; in the other words g*^^-' is comparable with the Coulomb 
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repulsion in the one-band Hubbard model. Hence, g^^^'/g^^'' ^ 10^. Then, by using Eqs. pT|) . ([T2|) and (|15l) we can 
calculate the ratio U/V. The result shows that U/V is proportional to 10 (g^^Vf'^') • Due to this reason U/V can 
be considerably larger than g^'^^ /g^^^- 

From the mathematical point of view, the value of U is not as much important as the value of the mean-field 
potential {Umf = It is easy to see that, in contrast to the BCS pairing potential, Umf is A-dependent (the 

fundamental feature of the presented model). So, the potential Umf depends on the temperature, V and U . If we 
set V and U the mean-field potential reaches the maximum value for T = K: U^^f = I^a/fIt-^q ~ where 
the symbol lA*^"-*! denotes the amplitude of the anomalous thermal average for T = K. In the case of the analyzed 



superconductors, we have 



(0) 
MF 



3t (see also subsection C). 



By using the Hamiltonian (fT4|) we calculate the anomalous Green function: 



((cktl 



)) = - 



V- 



^1 

6 I 



A 



(16) 



where iJk 



\/4 + iv + 



A|2. The obtained propagator is more complicated than the BCS Green function; 



peculiarly, we draw the readers' attention to the very intricate structure of the energy gap. The fundamental equation 
can be found in the form: 



(17) 



In order to calculate the thermodynamic properties we transform the momentum summation over an energy integration 
in Eq. pT)) . We notice that, in the case of three-dimensional system, where the electron density of states near the 

Fermi energy is constant, the mean- field 4EE interaction can be neglected because the value of U^^^p is very small. 

The situation changes dramatically for the two-dimensional system, where can be even greater than V. Then: 



1 = Vtr 



dep (e) 



tanh 



E 



(18) 



where: E = ye 
is given by [36j-[4l| 



-I- A2 



, Atot = Vtot I A| and Vtot = ^ + ¥|Ap. In the case of the square lattice the density of states 



p{e) = &iln|: 



(19) 



where h = -0M687t-^ and 62 = 21.17796i. 

In Figs, inland |3] we show the numerical solutions of Eq. (|18p for increasing values of C/|,/^. Analysis of the 



presented results allows one to state that only for C/|/^ < U^'' (where U^' is some characteristic value) the gap 

equation has one solution. Above Uq ^ at Tc open the two new branches of the energy gap. 

We notice that in the framework of the obtained mean-field model the 4EE interaction does not influence on the 
value of Tc; so, the critical temperature can be calculated by using the expression [soj : 



(0) 



B^C 



abye 



(20) 



where 



1 



In^ (2a) 



InM^l- 



- 2 



(21) 



and a = jn ~ 1.13 (7 is the Euler constant). In this case the isotope coefficient is small and can be calculated 



from: a = ^ In 



Returning to the central line of the thought, it can be easily seen, that for both regions of vfip the values of A^ot 

strongly increase when vf^p increases. For U^^f ^ ^ evolution of the gap parameter with the temperature is 
sharply different from the classical BCS prediction (see Fig. 2]). In particular for < T < the superconducting 
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FIG. 3: The dependence of Atot on the temperature for U'^pp < C/i?' = 0.2343t. We assume V — It and loq — O.Sf. 
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FIG. 4: The dependence of Atot on the temperature for U^^^p > Uq\ We assume V = It and ujq = 0.3t. The vertical Hne 
indicates a position of the critical temperature. For T > Tc the solid line represents the higher branch, whereas the perforated 
line corresponds to the lower branch. 



gap is very weak temperature dependent; this anomalous behavior is frequently observed in the cuprates I42l|. The 
another important results are presented for Tc < T < T**, where T** denotes the highest value of the temperature 
for which the non-zero solution of the gap equation exists. In this case we have two branches. In order to find out, 
for which of these solutions the thermodynamic potential is lower, the numerical calculations have been made. The 
detailed analysis of this complicated issue is presented in the Appendix C. As an example, in Fig. [5] we show the 
dependence of the difference of the thermodynamic potential between the non-zero gap state and the normal state 
{Aflvu) on the temperature. The obtained result proves, that the physical solution represents the higher branch; 
whereas the lower branch corresponds to the unstable state. 

Next, the temperature T** should be interpreted on the experimental background. 

First, we notice that the complicated mathematical structure of the order parameter (in general, it is the complex 
function) imposes two conditions on the existence of the superconducting state: (i) the amplitude of the order 
parameter has to differ from zero and (ii) the superconducting state has to exhibit the long-range phase coherence. 

The essential pointer, how should be interpreted the temperature T**, is connected with the experiments based 
on the Nernst effect @, 43 1. Namely, the Nernst signal above Tc strongly suggests that the superconductivity 



vanishes at Tc because the long-range phase coherence is destroyed by the thermally created vortices. Additionally, 
the experimental data have shown that the amplitude of the order parameter extends into the " normal" state to the 
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FIG. 5: The dependence of AQvu on the temperature for U^j^p — 1.2370t. We assume V — It and luq = 0.3t. The vertical hne 
indicates the position of the critical temperature. The solid line for < T < Tc represents the result for the superconducting 
solution; the solid line for Tc <T < T** corresponds to the higher branch. The perforated line represents the thermodynamic 
potential for the lower branch. We notice that A^lvu for the lower branch has the jump at Tc- This behavior is connected 
with the inversion of the solutions of Eq. (|18p in the considered region viz., for Vi and Ui higher than V2 and U2 we have 
Atot{Vi,Ui) < Atot{V2,U2) ■ 



temperature T/vb (the Nernst temperature). We notice that Tne is considerably much lower than the pseudogap 
temperature (T*). On the basis of presented experimental facts and the obtained theoretical results we assume that 

Before the comparison between the experimental data and the theoretical results obtained in the framework of toy 
model, we supplement the analytical approach. First, we notice that for T = the integral equation (|18p reduces to 
the algebraic equation: 

A[°l^2uJoe-^, (22) 
where the symbol A^^^ denotes the gap parameter at zero temperature and 



A2 ( &2 



Vtoth 



(23) 



By using the equation and expression (I^Hl) one can easily calculate the ratio: i?i = §^3^. In the classical BCS 
theory Ri is the universal constant of the model and [^ilgcg = 3.53. In the case of the BCS van Hove scenario Ri 
depends on the model parameters, however slightly ([i?i]niax ^ 4). The results obtained in the framework of the toy 
model are shown in Fig. HI We see that Ri is always bigger than in the BCS van Hove scenario and for sufficiently 
big values of U^j^p the ratio Ri achieves the physically acceptable values. 

Next, we consider the low-value behavior of the energy gap (Atot/ksTc << 1). In this case, the equation ^TE\\ 
should be rewritten in the form: 

9T/ +°° r'^o 1 
The kernel of the Eq. (|24|) may be expanded in powers of Af ot ■ 



1 1 ^It 



AL Ljl + e^ (^2 ^,2^2- 



(25) 



Next, we assume ujq +00, since luq >> ksTc- By using the lengthy but straightforward calculation we can 
transform the right-hand side of Eq. (IMl) into the algebraic form: 

l = b^Vtot[pi{T)-p2{T)Al,], (26) 
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FIG. 6: The dependence of the ratio Ri on U^j^p for different values of V. We assume ljq = O.St. We notice that the ratio Ri 
is plotted for different maximal values of U^j^p, since this quantity is also strong dependent on V. 



where 



and 



In' 



^0 

2fcsT 



(27) 



P2 (T) 



1 



irknT 



026 



(28) 



The symbol ^ denotes the Riemann zeta function, i) is defined by: -d (z) = (2j + 1) ^ In (2j + !)• If we omit the 

terms higher as |A|^ in Eq. ((26|) . the expression for the anomalous thermal average takes the form: 



n (T) ± ^rf {T)+T^{T)n (T) r^-i (T), 



(29) 



where: 



n (T) = -Upi (T) - (T) ^ 



(30) 



r2 (T) = -C/l/2p2 (T), 



(31) 



and 



rs (T) ^ 2 



Pi (T^)V^-7- 

Ol 



(32) 



We would like to point out that the expression (|29p can be used only in the framework of the generalized mean-field 
model where J7 0; in the case of the BCS van Hove scenario the formula generates the indeterminate form. 



In Fig. [7] we illustrate the temperature dependence of the gap parameter close to Tc for different values of 



(0) 
MF- 

In particular, we have shown the results obtained by using Eqs. (IT8)) and ((29)) : for comparison we calculate also the 
gap parameter in the BCS van Hove scenario. One can see that the approximate formula very well reconstructs 

the exact numerical results both for the small and large values of U^^p-, we strongly notice that for U^^p 
only the unstable state. 



> Uq^ exists 
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FIG. 7: The dependence of Atot on the temperature close to the transition temperature for the smaU and large values of U^^^p. 
We assume V = It and uo — 0.3t. The empty circles correspond to the exact numerical calculations (Eq. (|18|) V The solid line 
represents the calculation of Atot by using the formula ((29]). The crosses are obtained in the BCS van Hove scenario. 



C. The toy model and the experimental results 

The knowledge of the experimental values of Tc, T** and other parameters {t and ujq) enables the calculation of 

V and u'^^j'p for the real materials. In particular, the value of V is obtained by using the expression (I^H)) . Next, 

the potential U^^p is determined from Eq. (|18l) . We notice that, in principle, the quantities V and u'fip can be 
determined for the underdoped, optimally doped or overdoped samples; one can even take under consideration the 
influence of the disorder on V and U^mf- '^^^^ possible, since the strength of the presented approach lies in fact 
that the physical state of the high-Tc materials in the superconductivity domain is well reproduced by the values of 
Tc and T** . 

In the subsection (for selected superconductors), on the basis of V and C/^J^, we calculate the dependence of the 
ratio i?i on the hole density or on the doping; we consider also the influence of the disorder on Ri. We compare the 
theoretical predictions with the experimental results (if the experimental data is known) . In one case we plot openly 
the dependence of the energy gap on the temperature, since the right experimental data has existed in the literature. 

In particular, we take under consideration following families of the high-Tc materials: YBa2Cu307_y (YBCO) 
- for selected values of the hole density [p - (holes/Cu)), YBa2Cu307_y (YBCOo) - the disorder induced by elec- 
tron irradiation, which results in the creation of point defects such as Cu and O vacancies in the Cu02 planes, 
YBa2 (Cui_xZnx)3 Oy-y (Zn-YBCO) - the in-plane disorder caused by zinc, Yi_xPrxBa2Cu307_y (Pr-YBCO) - the 
out-of-plane disorder caused by praseodymium, NdBa2(Cui_xNix)307_y (Ni-NdBCO) - for selected values of x and 
y. Next, we consider Bi2Sr2CaCu208+y (Bi2212) - for selected values of the hole density, and Bi2Sr2Ca2Cu30io+y 
(Bi2223) - for the optimally doped case (OP). Finally, Pr2_xCexCu04_y (PCCO) - the example of electron-doped 
cuprate. 

The experimental data and the obtained theoretical results are collected in the Tab. H] Mathematically, the range 
of the Nernst region can be characterized by the quantity: 5** = {T** — Tc) /Tc- On the basis of experiments, we 
conclude that the Nernst region strongly expands if the hole density decreases; this is clear to see for YBCO and 
Bi2212, where [5**]^^^ = 7.08 and [5**]^^^ = 1.18 respectively; also for underdoped Ni-NdBCO, where [5**]^^^ = 3. 
In the case of the disorder which is induced in YBCO, the value of 5** considerably increases if the electron irradiation 
is used or yttrium is substituted by praseodymium; the disorder caused by zinc weakly influences on the value of 5** . 
In the presented analysis we consider also the electron-doped cuprate system PCCO for which 5** increases if doping 
decreases (from optimally to underdoped region). In overdoped region of PCCO the value of 5** can not be determined, 
since Tc and T** are experimentally indistinguishable. 

Now, we consider the obtained values of V and U^j^p. We notice that for real materials both V and U^^^p are 
significant and the following principle is in effect: if 5** is smaller than ~ 0.6 we have V > C/}5f' opposite case 

V < u'j^i^p; the especially high values of uj^^p can be observed in the strongly underdoped regime. 

By using the input parameters presented in Tab. U we compare the calculated theoretical values of the ratio Ri 
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TABLE I: The quantities S** , V and U];,'p calculated by using the mean values of Tc and T** . For YBCO the hole density 
p is estimated from the relation presented in [3]; for Bi2212 from the empirical formula Tc (p) /Tc,max = 1 — 82.6 {p — 0.16)^ 
[45| . For the superconductors Ni-NdBCO and Bi2223 the values of t and loo are unknown. In this case, we take t and loo for 
YBCO and Bi2212 respectively. 

Material Type t (mcV) Ref. ujo (meV) Ref. Tc (K) T" (K) Ref. S" V (t) (7^°'^ (t) 



YBCO 




p = 


0.062 


250 


[32] 


75 


m 


18.6 


150.2 


[46] 


7.08 


0.838 


2.847 






p = 


0.079 






45 


128.3 


1.85 


1.148 


2.045 






p = 


0.107 










60.5 


91.8 




0.52 


1.297 


1.171 






p = 


0.116 










64.1 


84.9 




0.32 


1.330 


0.969 






p = 


0.120 










66.5 


87.4 




0.31 


1.352 


0.974 






p = 


0.138 










80.6 


104.4 




0.30 


1.475 


1.062 






p = 


0.150 










90 


105 




0.17 


1.554 


0.898 






p = 


0.176 










92 


107 


[43] 


0.16 


1.571 


0.903 


YBCOe 




p = 


0.098 


250 


m 


75 


\21\ 


57 


85 


m 


0.49 


1.265 


1.104 






p = 


0.098 










45.1 


83.1 


0.84 


1.150 


1.288 






p = 


0.098 










24.2 


75 




2.10 


0.915 


1.564 






p = 


0.157 










92.6 


103 




0.11 


1.576 


0.796 






p = 


0.157 










79.5 


97.1 




0.22 


1.465 


0.930 






p = 


0.157 










48.6 


82.5 




0.70 


1.184 


1.213 


Zn-YBCO 




X = 


0.000 


250 


[32] 


75 




90 


104 ±5 


[48] 


0.16 


1.554 


0.876 






X = 


0.005 






84 


96 ±5 


0.14 


1.504 


0.813 






X — 


0.010 










79 


87 ±5 




0.10 


1.461 


0.700 






X — 


0.020 










67 


75 ±5 




0.12 


1.356 


0.674 


Pr-YBCO 




X - 


= 0.0 


250 


[32] 


75 


\21] 


89.7 


104.8 ± 5 


[49] 


0.17 


1.552 


0.899 






X - 


= 0.1 










83.8 


99.9 ±5 


0.19 


1.502 


0.908 






X - 


= 0.2 










68.2 


95 ±5 




0.39 


1.367 


1.096 






X = 


= 0.3 










50.2 


84.8 ± 5 




0.69 


1.200 


1.226 






X = 


= 0.4 










40.7 


79.9 ± 5 




0.96 


1.104 


1.313 


Ni-NdBCO 


X 


=0.00,y=0.0 


250 


- 


75 


- 


95 


115 ±20" 


[50] 


0.21 


1.596 


1.021 




X 


=0.03,y=0.0 










59 ±4 


80 ± 20" 


0.36 


1.283 


0.966 




X 


=0.06,y=0.0 










45 ± 5 


65 ± 20" 




0.44 


1.148 


0.940 




X 


=0.00,y=0.2 










53 ± 3.5 


75 ± 5 




0.42 


1.227 


0.981 




X 


=0.03,y=0.2 










20 ±7 


80± 10 




3.00 


0.858 


1.737 


Bi2212 




P = 


0.087 


350 


[51J,[52], 


80 




50 


108.9 ±5 


[43], [56]. 


1.18 


1.117 


1.508 






p = 


0.118 




[53]. 




[19], [5^, 


77.9 


130.3 ± 10 




0.67 


1.347 


1.442 






P = 


0.160 








[55]. 


90.6 


125.4 ± 5 




0.38 


1.444 


1.184 






p = 


0.202 










76.9 


105.5 ±5 




0.37 


1.339 


1.050 






P = 


0.219 










64.9 


85.5 ±10 




0.32 


1.243 


0.888 


Bi2223 




OP 


350 




80 




109 


135 


[43] 


0.24 


1.580 


1.077 


PCCO 




X = 


0.13 


380 


[33], [34], 


33 


[59], [60]. 


12 ± 1 


18 ± 1 




0.50 


0.795 


0.715 






X = 


-- 0.14 




157],[58]. 




19.5 ±1 


23 ± 1 


0.18 


0.947 


0.529 






X = 


0.15 










20 ± 1 


22 ± 1 






0.956 


0.174'' 






X = 


-- 0.17 










13 ± 1 


13.5 ± 1 






0.817 


0.002' 






X = 


= 0.19 










6.5 ± 1 


7± 1 






0.654 


~o' 



"The high value of the experimental error. 
''The value found based on IRi]^^^. 



with the experimental data. 

In Fig. [8]we show the ratio Ri as a function of the hole density for YBCO. The numerical results are denoted by the 
solid line with the open circles; the experimental data by the black filled symbols (see also the Tab. |lT]in Appendix 
D). It can be seen that, from slightly underdoped to overdoped crystals the presented model correctly predicts the 
values of Ri (taking under consideration the several approximations which have already been mentioned previously). 
For p 6 (0.07,0.135) the theoretical values of the ratio i?i are probably lower than the experimental data. However, 
the increasing of i?i with decreasing of the hole concentration is qualitatively correctly predicted. In the case of the 
strongly underdoped crystals {p < 0.07) the toy model suggests very high values of the ratio Ri which have to be 
experimentally checked. In the inset in Fig. [8] there is presented the open dependence of the critical temperature on 
the hole density with help of which the values of p for YBCO have been calculated 
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FIG. 8: The dependence of the ratio _Ri on p for YBCO. The sohd hne with the open circles represents the theoretical 
calculation. The black filled symbols correspond to experimental results obtained by: (a) - M. Sutherland et al. [s^], (b) - K. 
Nakayama et al. [H^Jc) - A. Kaminski et a/, [pi, (d) - M. Plate et al. 651, (e) - D.K. Morr et al. SjH.F. Fong et al. 
(f) - N.C. Yeh et al. (gll, (g) - V. Born et al. '6^(h) - H. Murakami et al. 70], (i) - H. Edwards et aL^, (j) - H. Edwards et 
al. [73 1 . The inset shows the dependence of the critical temperature on the hole density estimated from the relation presented 
in or 




0.000 0.005 0.010 0.015 0.020 

X 



FIG. 9: (A) The dependence of the ratio Ri on Tc for YBCOe; we consider the case p = 0.098 and p = 0.157. (B) The 
dependence of the ratio Ri on x for Pr-YBCO. (C) The dependence of the ratio Ri on x for Zn-YBCO. 



In Figs. [S](A)-(C) for YBCO we presented the influence of the disorder on the value of the ratio Ri. In Fig. [3] 
(A) , the dependence of Ri on Tc for the disorder induced by electron irradiation is shown. The two values of the hole 
concentration are considered: p = 0.098 and p — 0.157. In both cases the growing disorder, which is being manifested 
by the lowering values of Tc, causes very distinct increase of the ratio Ri. It is possible to observe the similar effect 
for the case, when the out-of-plane disorder is caused by praseodymium (Fig. [9](B)). However, the in-plane disorder 
which is induced by zinc in principle doesn't affect the value of the parameter Ri; see Fig. [5] (C). In Figs. [HI (A)-(C) 
we have shown only the theoretical predictions, since the experimental data not yet exist. 

The dependence of the ratio i?i on the hole density for Bi2212 is presented in Fig.fTU] Similarly as for YBCO, 
the theoretical results are denoted by the solid line with the open circles and the experimental data by the filled and 
half- filled symbols (the Tab. IIVI in Appendix D). Additionally, the region between the dotted lines represents the 
possible experimental values obtained by using the empirical formula [73|: 



Ri {p) = (15 ± 1) - (38±5)p. 



(33) 
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FIG. 10: The dependence of the ratio Ri on p for Bi2212. The sohd line with the open circles represents the theoretical 
calculation. The filled and half-filled symbols correspond to experimental results obtained by: (a) - Ch. Renner, et al. (b) 

- A. Hoffmann^t al. @, (c) - Y.G. Ponomarev, et al. [zl, (d) - T. Oki, et a/JTl, (e) - V.M. Krasnov, et al. [z3l, (f) - A.K. 
Gupta, et al. [tII, (g) - A. Kanigel, et al. 79], (h) - J.C. Campuzano, et al. [sfljT^K. Tanaka, et al. [sj, (i) - T. Nakaono, et 
al- m, (j) - M. Oda^ei al. [H, (k) - K. McElroy, et al. @, (1) - A. Matsuda, et al. [83], (m) - J.E. Hoffman, et al. (n) 

- C. Howald, et al. [13], (o) - H. Murakami, et al. [s^. The lines (p) are obtained by using the empirical relation (|33[l . The 
inset shows the dependence of the energy gap on the temperature; the filled region between the solid lines represents possible 
theoretical values of Atot', (r) - the experimental results [791 ]. 



We can notice that relation ([55)1 represents the linear least-squares fit to the experimental data and it is valid for 
0.12 < p < 0.24. On the basis of the presented results we conclude, that the toy model, in a wide range of the p 
values, very correctly reproduces the experimental data. It is possible to have reservation only for very low values of 
p, where some experimental data suggested extremely high values of the ratio The inset in Fig. [TU] presents the 
dependence of the energy gap on the temperature for Bi2212 with Tq = 67 K. In particular, the region between the 
solid lines corresponds to the possible theoretical values of Atot. Let us notice that theoretical results were set with 
the experimentally accuracy of T** (see also the Tab. In the inset we also show the experimental data obtained by 
A. Kanigel, et al. [79|]. It is easy to see that the agreement between theoretical and experimental results is excellent. 

However, for Bi2223 and Ni-NdBCO superconductors the values of t and luq are unknown and we take t and loq for 
Bi2212 and YBCO respectively, it is possible to receive good estimating of the real values of In the case of Bi2223 
we have Ri = 8.01. On the basis of the results presented in Appendix D (Tab. |V|) we see that our theoretical value 
is very close to the experimental data. For Ni-NdBCO the situation is more complicated. If we consider Ni-NdBCO 
superconductor with y=0 the experimental error for T** is too big, so it is impossible exactly determine Ri. However, 
if we take the mean values of Tc and T** the theoretical results indicate that the ratio Ri has the values 7.97, 10.20, 
11.78 respectively for x=0, x=0.03 and x=0.06. We notice that the experimental value of Ri for x=0 is equal to 7.3 
(see the Tab. lIIIl in Appendix D). For the case y=0.2 the experimental data are much more accurate and we have i?i 
equal to 10.98^1;^^ and 32.90t 

10.68 fo"" x=0.03 respectively. On the basis of above results we see that the 

last value of the ratio i?i can be extraordinary big (this result should be checked experimentally) . 

The general phase diagram of the high- Tc superconductors presents the global symmetry between the hole- and 
electron-doped materials [831 ■ First of all, in both cases, the antiferromagnetic phase has the similar Neel temperature 
(however for electron-doped cuprates, the antiferromagnetic phase is broader). Secondly, the superconducting phase 
for the hole- and electron-doped materials appears closely to antiferromagnetic phase, with the similar value of the 
optimal doping (x ^ 0.16). The distinct symmetry of the phase diagram strongly supports the view that the hole- 
and electron-doped superconductors should have an identical pairing mechanism. For that reason, the analysis of i?i 
for the selected electron-doped superconductor in the framework of toy model is very important. 

In Fig. [TT]we show the dependence of the ratio Ri on doping for PCCO. The solid line with open circles represents 
the theoretical calculation obtained by using the input parameters Tc and T**; the dotted line with open circles 
represents the theoretical calculation obtained by using the input parameters Tc and the appropriately selected ?7|5f- 
The black filled symbols correspond to experimental data. Based on Fig. [Til it is easy to see that the theoretical 
results very well reconstruct the experimental values of Ri . 
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FIG. 11: The dependence of the ratio _Ri on x for PCCO. The sohd hne with the open circles represents the theoretical 
calculation based on Tc and T**; the dotted line with open circles represents the theoretical calculation based on Tc and the 
appropriately selected C/mf- Th e black filled symbols correspond to experimental results obtained by: (a) - A. Biswas, et al. 
[93], (b) - A. Zimmers, et al. '91], (c) - Y. Dagan, et al. (d) - C.C. Homes, et al. [933, (e) - P- Fournier, et al. HI]. 



IV. SUMMARY 



In the study the simple microscopic model of the superconducting state that induce at high temperature was 
presented. The starting point of our considerations was the assumed statement: that a proper description of the 
superconducting condensate in the cuprates is possible only when the pairing mechanism would inseparably link 
together the strong electron correlations and the crystal lattice vibrations. The Hamiltonian proposed in the paper 
seems to be the most simple among the operators that satisfy the above postulate. 

In the study we have shown, that for the large values of the electron-electron-phonon coupling, the energy gap at 
the vicinity of the superconducting state existence weakly depends on the temperature and it vanishes above Tc (at 
the temperature that was interpreted as Nernst temperature). The key test of our model was the determination of 
the dependence of the ratio i?i on doping for the selected superconductors. The obtained agreement between the 
theoretical and experimental data seems to be extraordinary, when taking into account the simplicity of the considered 
Hamiltonian and used approximations. 

In our opinion, the achieved results clearly suggest that presented pairing mechanism should be seriously taken into 
consideration in the further researches which may lead to the most precise understanding of the properties of the high 
temperature superconducting state. 
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Appendix A: The exact expression for the self-energy matrix 

In the framework of the Eliashberg formalism, the Nambu spinors are defined in the following way: 

^^^(^y^J,^i^{ci ch), (Al) 
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where and ^'j^ satisfy the anticommutation rules: 

*k,<,]^ = V^o. (A2) 
In the presented notation, the terms of the Hamihonian ([T]) can be rewritten as follows: 

H^°^ = J2 ek*Lr3*k + J2 ^q^q^i' (A3) 

k q 

^^'^-Effk''^(q)*U^3*k0q, (A4) 

k 

and 

kk'ql kk'ql 



The matrix elements take the form: 



g(^'^)(q)^g(^)(q)+g(^)(q), (A6) 



where: 

fff (q) ^ E -i'^k 0) = E ffi'lk ' 



and 

,(2) 



N I 9, , , ' , (q, —1) for X = a 

fl(2f (q,l) = i '^^-l.k+l-qV^' y ^ 

^k,k ) I ^(2^ for x = 6. ^ ^ 



The bases of the matrices ta-td and tq — T3 (the Pauli matrices) are presented below: 

1 










TA = 


(J 









1 


TB = 














TC = 


(; 












TD = 




1 



2 (^0 + T3) , (A9) 

I (ri + iT2) , (AlO) 

I {ri - ^ra) , (All) 

I {ro - T3) . (A12) 

With the help of definition (jAip we have introduced the electron Green function Gk(jw„) = <^<^5'k|^'k^^ , where: 

ujn = (tt//?) (2n — 1) is the n-th Matsubara frequency and (3 = (ksT)"^ (kB is the Boltzmann constant). The electron 
propagator G\^{iuJn) can be written in the form: 

Gk(z..„) = ( t'^^l'V t'"'"'' 1 • (A13) 

V « C-kJcJ,^ »^-'„ << cLk4,|c_ki J 

We notice that the superconducting thermal average should be evaluated from the non-diagonal part of G]i{iujn), 
whereas the diagonal part determines the normal-state properties. It is essential, that the Hamiltonian (jlj in the 
matrix form enables to obtain the Dyson equation: 

Gk(ia;„) = Gok(i'^n) + Gok(i'^n)Mk(it^n)Gok(ia;„), (A14) 
where Gok(^Wn) denotes the unperturbed {g^^^ (q) — g'^^^, (q, 1) = 0) propagator: 

Gok(«t^«) = (iw„ro - ST^y^ , (A15) 

and the self-energy matrix has the form: 
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Mk(«w„) 



(A16) 



<< 



q q' 

q 



>>,: 



<< 



>>,: 



The equation (|A16P represents the formally exact expression for the self-energy matrix. The first-order contribution 
to Mk(zw„) comes from the thermal average () and the second-order contribution from the propagators in angular 
brackets. The complete analysis of the matrix self-energy is very difficult problem. However, one can demonstrate 
that, in the case of the lattice distortion's absence ((6q) = 0), the first-order contribution to the self-energy is equal 
to zero. In particular: 





(A17) 



where i/*-^' represents the fermionic parts in the Hamiltonian H'^'^K We notice that the non-diagonal part of the 
expression (|A17|) is proportional to the superconducting thermal average, whereas the diagonal part is connected with 
the average number of electrons per lattice site. 

The standard electron- phonon self-energy is represented by the second term in Eq. (|A16[) . which is proportional 

to ^ffk^'^^ (q)^ • On the basis of Eqs. (|A6p and (|A7|) it is easy to see that the value of the electron-phonon couphng 
can be directly increased by the electron-electron-phonon interaction. One should clearly emphasize the fact that the 
existence of the 17^^^ (q) element is the direct consequence of the form of the Hamiltonian ^ in the symmetric Nambu 
notation. From the physical point of view, the above result is highly nontrivial. On the other hand, the remaining 
terms also give the important contributions to the electron-phonon pairing and, on the Eliashberg equations level, 
they can not be neglected. 



Appendix B: The fold mean- field approximation of the 4EE Hamiltonian 

We have rewritten the interaction term in the Hamiltonian (1101) in the form: 



Hi, 



U 



UJO 



24iV3 

kk ql(7 



E 4-i<TCk<.\\q^ci:k^,_^c_k-<T, (Bi) 



where: = c' , c' , c , < c, / . By using the well known expression: AB c:^ (A) B + A (B) — (A) (B) , 

k Iqo- k +l+q-cr -k -l-qcT -k cr k -cr *= Jr' \ I ' \ I \ I \ n 

we have obtained: 

'^k'lqa - ^k'<TCk'+l+q_<,C^_k'-l-q£7 + ^k' +l+q<T^-k' (xV " ^kV^k' +l+qrT " (^2) 

The symbol Ako- is given by: Ako- = (c-kcrCk-o-)- In the next step, we have substituted (jB2[) into (|B1I) . The 
Hamiltonian can be rewritten as: 



Hi, 



TT ^° TT ^° 



^k'+l+q<TV-fT"-kl<Tf'-kV 24iV3 ^ ^kVV+l+q-£T"kl<T^_k'-l-q<T ^°'^> 



24iV3 

kk' qlfT kk' qlcr 



UJQ 



^ V A*, A . /i^'^ 

^k'+l+qcr^k'<T"-klo-' 



247V3 ^ k +l+q. 
kk qlCT 
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where: h\^y^ = c_j^_|_j_^cjj._[^CkcrC-k-cr- We notice that the linear terms with reference to the anomalous thermal 
average (Ako-) can exist only on this stage of the presented approximation. 

The operator (|B3p has still very complicated form (the first and second terms comprise the six fermion operators) . 
In order to solve the problem we have simplified the Hamiltonian (|B3I) again: 

hkl - ^-k<Tclk+i-a4-i<T + A* k+i^CkaC-k-<T - A_k^A-k+i^- (B4) 

Joining the expressions (IB3p and (jB4[) . we have obtained the sum of the terms, which are proportional to A^^, A^^ 
or A^^. Next, long in the form but straightforward calculations (in the presented mean- field scheme) give: 



Hi, 



U 
247V3 

U 
24/V3 

U 
127V3 



(B5) 



kk qla 



"^-k-lCTC-k-lo- + C_k-qo-cl^_qg. 



kqlcr 



kk qlo" 

On the basis of the expression (jB5[) . we see that in the obtained Hamlitonian exist only the terms, which are 
proportional to A^^ or A^^ (all operators with A^^ have reduced each other mutually) . We notice that the first term 
in Eq. (|B5I) can be neglected because it does not include the operators. In the case when we separate the momentums 
in the remained expression, the second term also can be neglected. The third term gives the interaction part of the 
Hamiltonian (jl3l) after using the mean-field approximation. 

It is clear that the used scheme of the simplification of the Hamiltoniam (jBll) is relatively simple. However, the 
presented approach lets us to analyse the physics of the considered system on the quantitative level. Thus, the 
agreement between obtained theoretical results and the experimental data will be an essential argument that the 
presented method is correct. The properties of the operator (jBl[) will be also studied diagramatically in the future. 
It will be done in order to improve eventually the theoretical description. 



Appendix C: The van Hove and generalized mean-field thermodynamic potential 



In the first step, we will calculate the thermodynamic potential in the framework of the BCS van Hove scenario. 
Next, we will generalize the results for the case U 0. 
For [7 = the thermodynamic potential is given by: 



n{v 



(-■) 



Pi Ha 



where the grand partition function has the form: 

z(f') =Tr 

The symbols Ha and Hb denote the following operators: 

ka 

and 



Hb 



kcr'^-k CT""k 



(CI) 



(C2) 



(C3) 



(C4) 



kk'o- 



The thermodynamic potential is readily found from the expression: 



d 
dV 



dV I3Z{V') 

+ 00 



Tr 



Integrate Eq. (|C5P from V = to y = y we obtain: 

1 



1. r 



dV (Hb) 



where Ay' =V A. The formula (|C6[) may be rewritten as follows: 

After substituting Eq. (for [/ = 0) into expression (IC7|) we find: 

Al^v = -f - 2 / rfep (e) (i; - e) 

+ -J^ dep{e)\n{l + e~^') 
A r^o 

- -J^ dep(e)ln(l + e-^^) 

The first integral in Eq. (jCSp is given by: 

Ji(Ay) = -2 / dep{e){E-s) 
3 

= hAl^Y^f^iAv), 
i=i 

where: 



1113 3 

2' 2' 2' 2' 2'^ \ A^ 



and 



/2 (Ay) ^ 



/3 (Ay) = In (62) 



- -ln(wo) 



^0 

Ay 



^0 
Ay 



Ay 



Ay 

Wo 



In 



- 1 



arcsinh 



_Wo_ 

Ay 
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FIG. 12: The dependence of AQv on the temperature. We assume V = It and lvq — 0.3t. The empty circles are obtained from 
Eq. (IC8|) with help of Eqs. (|C9|I and (|C14|I . Solid line represents the calculation of AQv using the formula (|C16|) . 



The symbol Fp_q (a; b; z) denotes the generalized hypergeometric function. 

The second integral in Eq. (|C8p describes the first temperature-dependent correction to the thermodynamic po- 
tential in the normal state: 



4 

h{T)^- d£p(e)ln(l + e-^-) 
P Jo 



(C13) 



The expression (|C13p can be rewritten by using the partial integration method. Since wo >> ksT, the obtained 
integral may be extended to infinity. In this way, we can find: 



/2(T) = -6i [k + In (62/3)] 



3/32 



(CM) 



+6iwo [In (a;o)-l]hi(l + e-'^"°; 



4 



The number k ~ 0.45403 is defined by: 

K = 1 + In (2) + 7 
'(971 (z) 



(C15) 



3 



dz 



fdjiiz) 
V dz 



where the symbol 7„ (z) is the generalized Stieltjes constant. 

We also notice that, by using the equations (|C7I) and (|26l) it is possible to derive explicit expression for Afly close 
to the transition temperature: 



Ani 



hi 



2P2 (T) 



biV 



The dependence of Afly on the temperature is shown in Fig. 1121 

The calculation of the thermodynamic potential for [/ 7^ is more difficult problem. In this case: 



where the grand partition function has the form: 

Z [V^u'j = Tr 
In Eq. (ICISP the Hamiltonian He is given by: 



P[ Ha- 



Hb — "24" 



(C16) 



(C17) 



(C18) 
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k' +l+q_aV-<TcLk'-l-qa'^-kVctk+l-crC-k-a- 



kk qlo- 

Now, we consider the total differential of the thermodynamic potential: 



, . dQ.(v\u'\ dn(v',U'\ 
(^'^)= 8V' '""^ dU' ' 



By using the method presented for the BCS van Hove scenario we can obtain: 



and 



-^^^ = -- {Hs) ^ -^|A (V,U) |^ 



(C19) 



(C20) 



(C21) 



(C22) 



Finally, the general evaluation of the thermodynamic potential requires the numerical analysis of the expression: 



Anvu = j^[n{v,u)-n{o,o)] 



(C23) 



7(00) l^r '^Jl''^^ +j^\^[v ,u)\Uu =- dx 



V\A {Vx, Ux) |2 + {Vx, Ux) I 
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Appendix D: Experimental values of Tc and the low-temperature superconducting gap 

In the Appendix we provide the list of the thermodynamic parameters vahies of high-Tc superconductors which 
have been obtained experimentahy. In particular, we have collected the data for Tc and the energy gap A[^1. We 
have also determined the doping level or stoichiometry of the materials and the values of Ri parameter. 



TABLE II: The experimental data for YBaaCusOr-y (YBCO). 



Type 


Tc (K) 


A[S (meV) 


Ri 


Ref. 


p = 0.079 
p = 0.111 
p = 0.166 
p = 0.184 


44 
62 
93.5 

89 


66 
71 
50 
37 


34.81 
26.58 
12.41 
9.65 


[62] 


p = 0.106 
p = 0.137 
p = 0.175 


60 
80 
92 


58 ±8.8 
45 ± 4.9 
34 ±3.1 


22.37 ±3.39 
13.04 ± 1.40 
8.55 ± 0.77 


[63j 


p = 0.099 
p = 0.110 
p = 0.140 
p = 0.160 


57.4 
61.8 
83 
93.2 


56.8 
54.2 
37.6 
33.8 


22.96 
20.36 
10.52 
8.42 


M 


p = 0.080 
p — u.uyo 
p = 0.160 
p = 0.184 


46.2 

J3D.O 

93.2 
88.9 


66.3 

/ 1 

49.9 
37.2 


33.27 

12.42 
9.7 


[65j 


p = 0.086 
p = 0.096 
p = 0.101 

p = O.i / i 


51.6 
56.3 
58 

no 


12.5 
14 
16.5 

on c; 
20.0 


5.62 
5.77 
6.6 

CC 1 o 

0.12 


[66], [67] 


p = 0.106 

p = 0.159 


60 
92.9 


~ 20 
18 


7.7 
4.5 


[68] 


p — 0.149 


89 


20 


5.2 




p = 0.153 
p = 0.153 


91 
91 


24-32 
~ 25" 


7.1 

6.4 


[70] 


p = 0.156 


92 


30 ±8 


7.6 


[71] 


p = 0.156 


92 


20 


5 


[72] 


"Tunneling 


on electrical field etched surface. 










TABLE III: The 


experimental data for NdBa2Cu307_y 


(NdBCO). 




Type 


Tc (K) 


^lol (meV) 


Ri 


Ref. 


y=o 


95 


30 


7.3 


[95] 
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TABLE IV: The experimental data for 612812 CaCu2 08+y (Bi2212) 



Type 


Tc (K) 


a'"' (meV) 

tot \ 1 




Ref. 


p = 0.125 


83 


44 ± 2 


12.30 ±0.56 


[42] 


p = 0.160 


92.2 


41.5 ± 2 


10.45 ± 0.5 




p = 0.208 


74.3 


34 ± 2 


10.62 ± 0.62 




p = 0.229 


56 


21 ± 2 


8.7 ± 0.83 




p — 0.131 


86 


32.2 ± 1.1 


8.7 ± 0.3 


f741 


p = 0.198 


81 


28.6 ± 1.4 


8.2 ± 0.4 




p = 0.189 


86 ± 4 


25 ± 1 


6.7 ± 0.5 


[75] 


p = 0.186 


87 ± 1 


26 ± 1 


8.2 ± 0.3 


[76] 


p = 0.165 


92 ± 3 


35 ± 1 


8.9 ± 0.7 


p = 0.134 


~ 87 


~ 45 


~ 12 




p = 0.127 


84 ± 4 


32 ± 0.5 


8.8 ± 0.6 




p = 0.123 


82 ±4 


33.5 ±0.5 


9.4 ±0.6 




p = 0.160 


92.5 


32.5 


8.15 


[771 


p = 0.181 


89 


25.8 


6.73 




p = 0.106 


70 


38 


12.6 


[781 


p = 0.191 


85 


30 ± 2 


8.2 ± 0.5 




p = 0.103 


67 


39.8 


13.79 


f79l 


p = 0.122 


80 


35.9 


10.42 


— u.uou 




fid ^ 


9Q 




^ — U.UOc' 


54.2 


61.2 


26.21 




p — U. -L -LU 


7^ 9 


47 S 


1 1 fi 




p — 0.115 


76.8 


50.1 


15.14 




p — 0.121 


ou.u 




1 '\ 97 




p = 0.133 


86.7 


43.5 


11.65 




p = 0.161 


92.2 


37.5 


9.43 




p = 0.186 


87 


31 


8.28 




p = 0.193 


83.8 


36.6 


10.14 




p = 0.201 


79.3 


25.8 


7.56 




— 205 


76.6 


34 


10.32 




p = 0.215 


69 


27.2 


9.13 




p = 0.100 


63 


40 


14.7 


[82] 


p = 0.130 


85 


33 


9.0 




p = 0.190 


85 


26 


7.1 




p = 0.100 


60 


36 ± 2 


13.9 




p — 0.140 


82 


34 ± 2 


9.6 


p = 0.160 


88 


32 ± 2 


8.4 




p = 0.210 


81 


27 ± 2 


7.7 




p = 0.110 


65 


62 


22.1 


f84l 


p = 0.130 


75 


48 ± 1 


14.9 




p = 0.150 


79 


43 ± 1 


12.6 




p = 0.180 


89 


36 ± 1 


9.4 




p = 0.190 


89 


33 ± 1 


8.6 




p = 0.110 


67 


55 ± 15 


19.1 


f85l 


p = 0.130 


85 


45 ± 12 


12.3 




p = 0.160 


89 


40 ± 10 


10.4 




p = 0.180 


89 


35 ±7 


9.1 




p = 0.220 


64 


22 ± 5 


8 




p = 0.120 


78 


50.2 


14.9 


[86] 


p = 0.160 


92 


43.7 


11 




p = 0.190 


85 


36.7 


10 




p = 0.120 


80 


42 ±2 


12.2 


[87] 


p = 0.120 


81 


40 


11.5 


[88] 
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TABLE V: The experimental data for Bi2Sr2Ca2Cu30io+y (Bi2223). 



Type 


Tc (K) 


(meV) 


Ri 


Ref. 


OP 


110 ±5 


36 ± 1.6 


7.6 ±0.5 


[75] 


UD 


109 


60 ±3 


12.8 


[96] 


OP 


111 


45 ±7 


9.4 





OP 109 ~ 37 ~ 7.88 [97] 



TABLE VI: The experimental data for Pr2-xCexCu04-y (PCCO). 



Type 


Tc (K) 


(mcV) 


Ri 


Ref. 


x=0.13 


12.2 


6.3 ±0.6 


11.99 ± 1.2 


[90] 


x=0.15 


21.6 


4.6 ±0.3 


4.95 ±0.36 


x=0.17 


11.8 


1.8 ±0.2 


3.46 ± 0.45 




x=0.15 
x=0.17 


21 
15 


(4.3-5.4) 
(3-3.6) 


(4.7-6) 
(4.7-5.6) 


[91] 


x=0.15 


19 ±1 


3.25 


4 ±0.4 


[92] 


x=0.16 


16 ± 1 


2.6 


3.8 ± 0.5 


x=0.17 


13 ± 1 


1.3 


2.3 ± 1.3 




x=0.18 


11 ± 1 


1.0 


2.1 ± 1.5 




x=0.19 


8 ±0.4 


0.9 


2.6 ±0.9 




x=0.15 


20 


4.3 


5 


[93] 


x=0.13 


10 


6.5 


15.09 


[94J 


x=0.15 


23 


4.4 


4.4 


x=0.17 


13 


1.75 


3.12 
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